Advances in Biomedical Engineering Research (ABER) Volume 1 Issue 3, September 2013 



www.seipub .org/ aber 



Heart Rate Variability in Essential 
Hypertension Patients with Different Stages 
by Nonlinear Analysis: A Preliminary Study 

Ping Shi, Hongliu Yu* 

School of Medical Instrument and Food Engineering, University of Shanghai for Science and Technology, Shanghai, 
200093, China 

*garendon@l 63 .com 



Abstract 

Little is known about the nonlinear dynamics of heart rate 
variability (HRV) in hypertension patients with different 
stages. The aim of this study was to investigate the difference 
of HRV in hypertension patients by nonlinear methods. The 
patients with stage 2 hypertension (n=10) and stage 3 
hypertension (n=12) were included in this study, in which 
the Allan factor (AF), approximate entropy (ApEn), Renyi 
entropy (REn), largest Lyapunov exponents (LLE) and 
wavelet-transform standard deviation (WTSD) as nonlinear 
characteristics of HRV were examined as well. The results 
showed that there are significant differences between stage 2 
hypertension patients and stage 3 hypertension patients in 
parameters of AF, ApEn, REn and LLE. The present findings 
indicated that nonlinear parameters are sensitive to 
discriminate hypertension patients with different stages, and 
patients with stage 3 hypertension have depressed nonlinear 
behavior compared with patients with stage 2 hypertension. 
Duplicated research with larger samples is needed in the 
future work. 
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Introduction 

Most population based studies confirm that 
hypertension profoundly increases an individual's risk 
of various cardiovascular consequences. Heart 
rate variability (HRV) has emerged as a practical, 
noninvasive tool to quantitatively investigate cardiac 
autonomic dysregulation in hypertension. HRV 
characterizes the variations in the intervals between 
consecutive heartbeats, that is, the variations in 
duration of R-R intervals between consecutive QRS 
complexes, measured by electrocardiograph (ECG). 
Assessment of HRV may provide quantitative 
information on the modulation of cardiac vagal and 
sympathetic nerve input. HRV analysis is a recognized 
tool for the estimation of cardiac autonomic 
modulations(Task Force 1996). Previous studies have 



found that reduced HRV is a powerful and 
independent predictor of an adverse prognosis in 
patients with cardiovascular disease (Bauer et al. 2008, 
Huikuri and Stein 2012). 

Numerous studies have assessed the relationship 
between hypertension and HRV(Garcia-Garcia et al. 
2012, Schroeder et al. 2003, Singh et al. 1998, 
Terathongkum and Pickler 2004, Yoo et al. 2011). 
Linear methods, including time- and frequency- 
domain analysis, are applied in these studies. Actually, 
nonlinear phenomena involved in the genesis of 
HRV are determined by complex interactions of 
haemodynamic, electrophysiological and humoral 
variables, as well as by autonomic and central nervous 
regulations. It has been speculated that analysis of 
HRV based on the methods of nonlinear dynamics 
might elicit valuable information for the physiological 
interpretation of HRV and for the assessment of the 
risk of sudden death (Task Force 1996). Nonlinear 
dynamical HRV analysis, based on nonlinear 
mathematics and chaos theories, has been made so as 
to improve the understanding of the underlying 
physiological processes of the autonomous nervous 
system (Hoyer et al. 1997, Kim et al. 2005, Krstacic et al. 
2012, Pivatelli et al. 2012). These methods do not 
attempt to assess the actual magnitude of HRV but to 
describe the complexity or fractal dynamics of the R-R 
interval time series. These methods may provide 
information beyond the conventional time and 
frequency domain methods (Makikallio et al. 1996, 
Voss et al. 1996) . 

In the present study, the nonlinear components of 
HRV were computed by using several methods, 
including Allan factor (AF), approximate entropy 
(ApEn), Renyi entropy (REn), largest Lyapunov 
exponents (LLE) and wavelet-transform standard 
deviation (WTSD). These methods differ from 
traditional methods, because they typically detect 
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qualitative rather than quantitative properties from HR 
time series. Also, these measures have been recognized 
by most researchers to examine the nonlinear 
characteristics of HRV (Ashkenazy et al. 1998, Beckers 
F. et al. 2001, Hagerman et al. 1996, Lake 2006, Turcott 
and Teich 1996). Most studies indicated that decreased 
nonlinear behaviour of heart rate is associated with 
pathological states. 

To our knowledge, few studies have examined the 
nonlinear characteristic of HRV in patients with 
hypertension, and none investigated the differences 
between patients who were in different hypertension 
stages. The aim of this study was to investigate the 
HRV difference between essential hypertension patients 
with stage 2 hypertension and stage 3 hypertension by 
using nonlinear components of HRV. 

Material and methods 

Subjects 

Essential hypertension patients were recruited in this 
study. 10 patients with stage 2 hypertension (SBP/DBP: 
160-179/100-109) and 12 patients with stage 3 
hypertension (SBP/DBP:>180/110) were investigated. 
The information of the subjects is presented in Table 1. 
All subjects participating in this study gave written 
informed consent. The local Ethical Committee 
approved the protocol of the study. 

TABLE 1 THE INFORMATION OF THE SUBJECTS. 





Stage 2 hypertension 


Stage 3 hypertension 


M/F 


2/8 


8/4 


Age (year) 


62.6±9.3 yr 


68.6±7.3 yr 


Height (cm) 


162±5 


166±4 


SBP(mmHg) 


171±10 


193±7 


DBP(mmHg) 


105±8 


125±9 



Study Protocol 

The ECG signals were collected while the patients 
were in a sitting position after resting for at least 20 
minutes. The measurements were taken in the morning 
and at the same room by one trained research assistant 
according to a standardized method. 8 min ECG 
recordings were continuously collected (1,000 samples 
per second) using a Powerlab Acquisition System 
(ADInstruments Corporation, Sydney, Australia). Beat- 
to-beat cardiac interval values, i.e. heartbeat series, 
were automatically measured for each sinus beat and 
exported for further analysis. 

There are many factors that may influence the HRV, 
including circadian rhythm, body position, activity 
level prior to recording, medication and breathing 
condition (Task Force 1996). For that reason, special 



A(7> 



precautions were taken to maintain similar condition, 
such performing both recordings in the morning and 
in a sitting body position after an adaptation time of at 
least 15 minutes. 

Nonlinear analysis of HRV 

1) Allan Factor 

The Allan factor is the ratio of the event-number 
Allan variance to twice the mean: 

2E{NjT)} (1) 
The Allan variance, as opposed to the ordinary 
variance, is defined in terms of the variability of 
successive counts (Thurner et al. 1997, Turcott and 
Teich 1996). 

2) Approximate Entropy 

The approximate entropy, a measure of system 
complexity introduced by Pincus (Pincus S.M. 1991), 
has been widely applied in cardiovascular studies 
(Ho et al. 1997, Kim et al. 2005, Lin et al. 2001, 
Perkiomaki et al. 2002). Two parameters m and r 
must be chosen prior to the computation of 
approximate entropy, where m specifies the pattern 
length, and r is the effective filter. Here also, one 
has to compute the correlation integral C m (r) 
(with embedding dimension m and time lag 1). This 
measure is finally obtained as follows 

1 



ApEn{m,r,L)- 



L-m 
1 



X^logC" +1 (r). 

I,tr 1 iogc;»(r) 



(2) 



L-m+1 

Thus approximate entropy quantifies the (log) 
likelihood that sets of patterns that are close remain 
close on next incremental comparison. Smaller 
values of approximate entropy imply a greater 
likelihood that certain patterns of measurements 
will be followed by similar measurements. If the 
time series is highly irregular, the occurrence of 
similar patterns in the future is less likely. Higher 
values of approximate entropy indicate a more 
complex structure in the time series. For this study, 
m is set to be 2 and r is set to be 15% of the standard 
deviation of each time series. These values are 
selected on the basis of previous studies indicating 
good statistical validity for approximate entropy 
within these variable ranges (Pincus S. M. and 
Goldberger 1994). 

3) Renyi Entropy 

The Renyi entropy of order a , where a > and 
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a ^ , is defined as 



I- a 



(3) 



Here, X is a discrete random variable with possible 
outcome 1, 2,..., n and corresponding probabilities 
/?, = Pr(X = i) for i=l,...,n, and the logarithm is 

base 2. If the probabilities are Pi = 1/ ra for all 
i=l,...,n, then all the Renyi entropies of the 
distribution are equal: H a (X ) = log n . m general, 
for all discrete random variables X, H a (X) i s a 
non-increasing function in a ■ 

Applications often exploit the following relation 
between the Renyi entropy and the p-norm: 

H a(X)- T ^lOg(\\X\\ a )_ (4) 

Here, the discrete probability distribution X is 
interpreted as a vector in 9?" with X. — p l >0 



and 



X ; =1 



4) Largest Lyapunov Exponents 

Lyapunov exponents represent the average 
exponential separation of initially nearby points. A 
positive Lyapunov exponent indicates that the 
initially neighboring trajectories diverge 
exponentially, which means that the system is 
sensitively dependent on initial conditions, the 
distinctive property of a chaotic system. An 
m-dimensional system has m Lyapunov exponents. 
Among those Lyapunov exponents, the largest 
Lyapunov exponent is enough to characterize the 
chaotic system in physiological applications 
(Radhakrishna et al. 2000). Our calculation is based 
on the algorithm developed by Rosenstein et al 
(Rosenstein et al. 1993), which is appropriate for 
small data sets. 

5) Wavelet-Transform Standard Deviation 

A dyadic discrete wavelet transform for the 
heartbeat sequence {r, } may be defined as: 



(5) 



1 L ~' 

W m,n M = ~j= Z WO'/™ " n) 

A simple measure that can be derived from the 
wavelet transformation is the standard deviation of 
the wavelet transform as a function of scale: 



<r wm M = [e[w^ (m) - E[w^ n (m)] 2 



(6) 



where the expectation is taken over the process of 
RR intervals, and is independent of n. It is readily 
shown that E[W mn (m)] = for all values of 
so that a simplified form for the 



m 



wavelet-transform standard deviation emerges: 



cr w Jm) = {E[W m Jm) 2 ]} 



(7) 



Statistical Analysis 



All available data were included in statistical analyses 
performed using the SPSS version 11 (Chicago, IL, 
USA). All data were presented as mean ± standard 
deviation (SD) and 95% confidence interval (CI) for the 
mean. The significance of difference between groups 
was compared using one-way ANOVA followed by 
post-hoc analysis. P less than .05 was considered 
statistically significant. 

Results 

Table 2 shows the comparison of nonlinear parameters 
between patients with stage 2 hypertension and stage 3 
hypertension. The values of mean±SD and 95% CI for 
the mean in nonlinear parameters are presented. 
Compared to patients with stage 2 hypertension, 
patients with stage 3 hypertension have smaller value 
of mean±SD and lower 95% CI for the mean for all 
nonlinear parameters. Specifically, in stage 3 
hypertension patients, parameters of AF significantly 
decreased (P=0.010), as well as ApEn (P=0.008), REn 
(P=0.006), and LLE (P=0.013). No significant differences 
were found for parameters of WTSD between patients 
with stage 2 hypertension and stage 3 hypertension 
(P>0.05). 



TABLE 2 COMPARISON BETWEEN PATIENTS WITH STAGE 2 HYPERTENSION AND STAGE 3 HYPERTENSION. VALUES ARE 95% CI FOR THE MEAN 





Stage 2 hypertension 
(n=10) 


Stage 3 hypertension 
(n=12) 


P 




Meanistandard deviation 


95% CI 


meanistandard deviation 


95% CI 




AF 


0.039±0.062 


0.027- 0.051 


0.022±0.013 


0.013- 0.029 


0.010 


ApEn 


1.64±0.62 


1.071- 2.21 


0.803±0.324 


0.63 -1.13 


0.008 


REn 


2.11+0.563 


1.70-2.51 


1.39±0.43 


1.13- 1.72 


0.006 


LLE 


0.134±0.042 


0.09-0.15 


0.082±0.039 


0.05- 0.10 


0.013 


WTSD 


32.6±19.7 


16.4 -48.9 


26.6±17.1 


9.8 - 38.5 


ns 



AF: Allan factor; ApEn: Approximate entropy; REn: Renyi entropy; LLE: Largest Lyapunov Exponents; WTSD: Wavelet-Transform Standard 
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FIG. 1 DISTRIBUTION OF VALUES FOR NONLINEAR PARAMETERS IN PATIENTS WITH STAGE 2 AND STAGE 3 HYPERTENSION. THE 

SIGNIFICANT DIFFERENCES (P<0.05) WERE MARKED WITH "*" 



In Fig. 1, the values of nonlinear parameters are 
presented, and box-and-whisker plots are applied to 
present the median and the values from the lower to 
upper quartile (25 to 75 percentile). Patients with stage 
3 hypertension have lower values in parameters of AF, 
ApEn, REn and LLE (all P<0.05) compared to patients 
with stage 2 hypertension. No statistically significant 
differences were found in parameters of WTSD 
between groups. 

Discussion 

Fluctuations in the heartbeats series represent the 
heart's ability to respond to naturally occurring 
physiologic variations in cardiovascular activity and 
regulatory impulses from various cardiovascular 
control systems. When the research focus is on 
dynamic changes of autonomic tone andor prediction 
of subsequent physiologic outcome, HRV is preferred 
over heart rate as an index of autonomic cardiac 
response. It is recognized that heartbeats fluctuations 
stem from dynamics of chaotic nature and exhibit 
short range correlation governed by deterministic laws. 
In addition, the predictive value including the time 
and frequency domain measures of HRV, is relatively 
low in some studies for risk stratification (Odemuyiwa 
et al. 1991, Redwood et al. 1997). In order to get a 
better understanding of heart rate regulation, methods 
based on nonlinear mathematics and chaos theories for 
risk stratification in clinical practice are necessary and 



required. Altered heart rate dynamics has been 
demonstrated in patients with chronic heart failure 
(Ho et al. 1997) myocardial infarction (Makikallio et al. 
1996) and a propensity for ventricular arrhythmias 
(Makikallio et al. 1999). No studies investigate the 
changes of heart rate dynamics in Essential 
Hypertension Patients. 

AF is a useful measure of HRV (Turcott and Teich 
1996). An excess above unity reveals that a sequence is 
less ordered than a homogeneous Poisson point 
process, while values below unity signify sequences 
which are more ordered. The lower AF in patients with 
stage 3 hypertension demonstrated the higher ordered 
heartbeats series, indicating the reduced HRV. 

ApEn and REn quantify the entropy of the system 
(Pincus S.M. 1991, Lake 2006). Entropy measures the 
degree to which the occurrence of a value depends on 
its predecessors in the input. The more regular and 
predictable the heartbeats series is, the lower the value 
of the entropy will be. On the other hand, a more 
random heartbeats series will lead to higher entropy 
values. Generally, healthy person have higher values 
of entropy compared to patients with impaired 
cardiovascular system(Makikallio et al. 2000). ApEn is 
an index of overall complexity and predictability of 
heart beat series. REn, which is the measures of heart 
rate gaussianity, was utilized to enhance and 
complement the analysis of biological signals using 
ApEn (Lake 2006). The Gaussianity of heartbeats series 
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is a measure that can be linked to physiological 
complexity. Lower ApEn and REn were observed in 
patients with stage 3 hypertension compared with 
patients with stage 2 hypertension, indicating that the 
more severe hypertension lowered the complexity of 
heartbeats series. Previous studies showed that heart 
diseases presented changes in the complexity of 
heartbeats series, assessed by ApEn. Vikman et al. 
(Vikman et al. 1999) revealed reduced ApEn preceding 
spontaneous onset of paroxysmal atrial fibrillation in 
patients without structural heart disease. ApEn is 
mostly linked with vagal modulation (Beckers F. et al. 
2001). The present findings suggested the reduced 
vagal nervous activity in patients with stage 3 
hypertension. 

The LLE quantifies the sensitivity of the system to 
initial conditions. This is characterized by the average 
rate of divergence of two neighboring trajectories in 
phase space and gives a measure of predictability. A 
negative value implies that the orbits approach a 
common fixed point. A positive value can be taken as a 
definition of chaos provided that the system is known 
to be deterministic(Rosenstein et al. 1993). Larger 
values of the LLE indicate more complex behavior. The 
smaller LLE in patients with stage 3 hypertension 
indicated the reduced HRV compared to the patients 
with stage 2 hypertension, suggesting the more severe 
impairment of autonomic nervous system in patients 
with stage 3 hypertension. However, in previous 
studies, this index was reported to be uncorrelated to 
autonomic modulation (Hagerman et al. 1996) . 

WTSD has recently been shown to be quite valuable 
for HRV analysis (Ashkenazy et al. 1998, Thuner et al. 
1998). The depression of WTSD is likely associated 
with the impairment of autonomic nervous system 
function. Although no significant differences were 
found between two groups in the present study, the 
patients with stage 3 hypertension did have a lower 
value of WTSD, indicating the reduced baroreflex 
sensitivity. 

Limitations 

Despite promising results in this study, there are some 
limitations to consider. First, this study was performed 
with small populations. Even nonlinear heart rate 
dynamic seem to be potential tools for risk 
stratification in hypertension patients, replication 
research with larger samples is needed to confirm the 
present findings. Second, several factors limit the 
clinical utility of nonlinear HRV analysis. Physiological 



factors, such as sex (Goldberger et al. 2002) and ageing 
(Beckers F. et al. 2006), have been shown to influence 
nonlinear heart rate dynamics. Although the two 
groups did not differ significantly in age and body 
weight, the difference of gender distribution in this 
preliminary study prompted the future investigation 
regarding gender difference when altered heart rate 
dynamics is interpreted. 

Conclusion 

In recent years, new dynamic methods of HRV 
quantification have been used to uncover nonlinear 
fluctuations in heart rate that are otherwise not 
apparent. The nonlinear variations would enable the 
cardiovascular system to respond more quickly to 
changing conditions. This study has investigated the 
nonlinear characteristics of HRV in patients with stage 
2 and stage 3 hypertension. Patients with stage 3 
hypertension have depressed nonlinear behavior of 
heartbeat series compared with those with stage 2 
hypertension. This can be related to the general 
concept of decreasing autonomic modulation due to 
high blood pressure. 

Although methods derived from nonlinear dynamics 
have been successfully applied to many scientific 
disciplines, the physiological background underlying 
nonlinear heart rate dynamics is far from being fully 
understood. Therefore, better understanding of the 
physiology of altered heart rate dynamics is required 
before these methods will gain any significant position 
in clinical practice. 

ACKNOWLEDGEMENT 

We thank Prof. Yi-Sheng Zhu, from Shanghai Jiao 
Tong University, Shanghai, China, for his expert 
technical assistance. 

REFERENCES 

Ashkenazy Y, Lewkowicz M, Levitan J, Moelgaard H, Bloch 
Thomsen PE, Saermark K. " Discrimination of the healthy 
and sick cardiac autonomic nervous system by a new 
wavelet analysis of heartbeat intervals. " Fractals 6 (1998): 
197-203. 

Bauer A, et al. "Heart rate turbulence: standards of 
measurement, physiological interpretation, and clinical 
use: International Society for Holter and Noninvasive 
Electrophysiology Consensus. " J Am Coll Cardiol 52 
(2008): 1353-1365. 



37 



www. seipub .org/aber 



Advances in Biomedical Engineering Research (ABER) Volume 1 Issue 3, September 2013 



Beckers F, Ramaekers D, Aubert AE. "Approximate entropy 
of heart rate variability: validation of methods and 
application in heart failure. " Cardiovasc Eng 1 (2001): 
177-182. 

Beckers F, Verheyden B, Aubert AE. Aging and nonlinear 
heart rate control in a healthy population. American 
Journal of Physiology-Heart and Circulatory Physiology 
290(2006): H2560-H2570. 

Garcia-Garcia A, Gomez-Marcos MA, Recio-Rodriguez JI, 
Patino-Alonso MC, Rodriguez-Sanchez E, Agudo-Conde 
C, Garcia-Ortiz L. "Office and 24-hour heart rate and 
target organ damage in hypertensive patients. " BMC 
Cardiovasc Disord 12 (2012): 19. 

Goldberger AL, Amaral LAN, Hausdorff JM, Ivanov PC, 
Peng CK, Stanley HE. Fractal dynamics in physiology: 
Alterations with disease and aging. Proceedings of the 
National Academy of Sciences of the United States of 
America 99(2002): 2466-2472. 

Hagerman I, Berglund M, Lorin M, Nowak J, Sylven C. 
"Chaos-related deterministic regulation of heart rate 
variability in time- and frequency domains: effects of 
autonomic blockade and exercise. "Cardiovasc Res 31 
(1996): 410-418. 

Ho KK, Moody GB, Peng CK, Mietus JE, Larson MG, Levy D, 
Goldberger AL. "Predicting survival in heart failure case 
and control subjects by use of fully automated methods 
for deriving nonlinear and conventional indices of heart 
rate dynamics. " Circulation 96 (1997): 842-848. 

Hoyer D, Schmidt K, Bauer R, Zwiener U, Kohler M, Luthke 
B, Eiselt M. "Nonlinear analysis of heart rate and 
respiratory dynamics. IEEE Eng Med Biol Mag 16 (1997): 
31-39. 

Huikuri HV, Stein PK. "Clinical application of heart rate 
variability after acute myocardial infarction. "Front 
Physiol 3 (2012): 41. 

Kim WS, Yoon YZ, Bae JH, Soh KS. "Nonlinear characteristics 
of heart rate time series: influence of three recumbent 
positions in patients with mild or severe coronary artery 
disease. "Physiol Meas 26 (2005): 517-529. 

Krstacic G, Parati G, Gamberger D, Castiglioni P, Krstacic A, 
Steiner R. "Heart rate variability and nonlinear dynamic 
analysis in patients with stress-induced 
cardiomyopathy. "Med Biol Eng Comput 50(2012): 



1037-1046. 

Lake DE. "Renyi entropy measures of heart rate 
Gaussianity. " IEEE Trans Biomed Eng 53 (2006): 21-27. 

Lin LY, Lin JL, Du CC, Lai LP, Tseng YZ, Huang SK. 
"Reversal of deteriorated fractal behavior of heart rate 
variability by beta-blocker therapy in patients with 
advanced congestive heart failure." J Cardiovasc 
Electrophysiol 12 (2001): 26-32. 

Makikallio TH, Tulppo MP, Seppanen T, Huikuri HV. 
"Analysis of nonlinear heart rate dynamics in cardiac 
arrhythmias." Herzschr Elektrophys 11 (2000): 131-138. 

Makikallio TH, Seppanen T, Niemela M, Airaksinen KE, 
Tulppo M, Huikuri HV. Abnormalities in beat to beat 
complexity of heart rate dynamics in patients with a 
previous myocardial infarction. J Am Coll Cardiol 28 
(1996) : 1005-1011. 

Makikallio TH, Hoiber S, Kober L, Torp-Pedersen C, Peng CK, 
Goldberger AL, Huikuri HV. Fractal analysis of heart rate 
dynamics as a predictor of mortality in patients with 
depressed left ventricular function after acute myocardial 
infarction. TRACE Investigators. TRAndolapril Cardiac 
Evaluation. Am J Cardiol 83(1999): 836-839. 

Odemuyiwa O, Malik M, Farrell T, Bashir Y, Poloniecki J, 
Camm J. Comparison of the predictive characteristics of 
heart rate variability index and left ventricular ejection 
fraction for all-cause mortality, arrhythmic events and 
sudden death after acute myocardial infarction. Am J 
Cardiol 68(1991): 434-439. 

Perkiomaki JS, Zareba W, Badilini F, Moss AJ. "Influence of 
atropine on fractal and complexity measures of heart rate 
variability. " Ann Noninvasive Electrocardiol 7(2002): 
326-331. 

Pincus SM. "Approximate entropy as a measure of system 

complexity." In Proceedings of the National Academy of 

Sciences of the U.S.A. (1991): 2297-2301. 
Pincus SM, Goldberger AL. "Physiological time-series 

analysis: what does regularity quantify? "Am J Physiol 

266 (1994): H1643-1656. 
Pivatelli FC, et al. "Sensitivity, specificity and predictive 

values of linear and nonlinear indices of heart rate 

variability in stable angina patients. "Int Arch Med 5 

(2012): 31. 

Radhakrishna RK, Dutt DN, Yeragani VK. "Nonlinear 



38 



Advances in Biomedical Engineering Research (ABER) Volume 1 Issue 3, September 2013 



www.seipub .org/ aber 



measures of heart rate time series: influence of posture 
and controlled breathing. "Auton Neurosci 83 (2000): 
148-158. 

Redwood SR, Odemuyiwa O, Hnatkova K, Staunton A, 
Poloniecki I, Camm AJ, Malik M. Selection of dichotomy 
limits for multifactorial prediction of arrhythmic events 
and mortality in survivors of acute myocardial infarction. 
Eur Heart J 18(1997): 1278-1287. 

Rosenstein MT, Collinsa JJ, De Lucaa CJ. "A practical method 
for calculating largest Lyapunov exponents from small 
data sets. "Physica D 65 (1993): 117-134. 

Schroeder EB, Liao D, Chambless LE, Prineas RJ, Evans GW, 
Heiss G. "Hypertension, blood pressure, and heart rate 
variability: the Atherosclerosis Risk in Communities 
(ARIC) study. " Hypertension 42 (2003): 1106-1111. 

Singh JP, Larson MG, Tsuji H, Evans JC, O'Donnell CJ, Levy 
D. "Reduced heart rate variability and new-onset 
hypertension: insights into pathogenesis of hypertension: 
the Framingham Heart Study. " Hypertension 32 (1998): 
293-297. 

Task Force of the European Society of Cardiology and the 
North American Society of Pacing and Electrophysiology. 
"Heart rate variability: standards of measurement, 
physiological interpretation and clinical use." Circulation 
93 (1996): 1043-1065. 

Terathongkum S, Pickler RH. "Relationships among heart 
rate variability, hypertension, and relaxation techniques." 



J Vase Nurs 22 (2004): 78-82; quiz 83-74. 

Thuner S, Feurstein MC, Teich MC. "Multiresolution wavelet 
analysis of heartbeat intervals discriminates healthy 
patients from those with cardiac pathology. "Phys Rev 
Lett 80 (1998): 1544-1547. 

Thurner S, Lowen SB, Feurstein MC, Heneghan C, 
Feichtinger HG, Teich MC. "Analysis, synthesis, and 
estimation of fractal-rate stochastic point processes." 
Fractals 5 (1997): 565-595. 

Turcott RG, Teich MC. "Fractal character of the 
electrocardiogram: distinguishing heart-failure and 
normal patients. "Ann Biomed Eng 24 (1996): 269-293. 

Vikman S, Makikallio TH, Yli-Mayry S, Pikkujamsa S, 
Koivisto AM, Reinikainen P, Airaksinen KE, Huikuri HV. 
"Altered complexity and correlation properties of R-R 
interval dynamics before the spontaneous onset of 
paroxysmal atrial fibrillation. "Circulation 100 (1999): 
2079-2084. 

Voss A, Kurths J, Kleiner HJ, Witt A, Wessel N, Saparin P, 
Osterziel KJ, Schurath R, Dietz R. The application of 
methods of non-linear dynamics for the improved and 
predictive recognition of patients threatened by sudden 
cardiac death. Cardiovasc Res 31(1996): 419-433. 

Yoo CS, Lee K, Yi SH, Kim JS, Kim HC. "Association of heart 
rate variability with the framingham risk score in healthy 
adults. " Korean J Fam Med 32 (2011): 334-340. 



39 



